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Abstract 

About 3fT deviation from the standard model prediction of muon anomalous mag- 
■ netic moment (muon g-2) has been reported. We consider new physics beyond the 

standard model which has new Yukawa interactions with muon. We compute new 
contributions to muon g-2 and corrections to electroweak observables, and show the 
consistent region of parameter space. We find that in a simple model where the chiral- 
ity flip of muon occurs only in the external muon line in one-loop muon g-2 diagrams, it 
is necessary to introduce the relatively large new Yukawa coupling and the electroweak 
scale new particles. On the other hand, in a model where the chirality flip can occur in 
the internal fermion line of one loop muon g-2 diagrams, we can obtain favorable g-2 
contributions without large Yukawa coupling, and they are consistent with the precision 
electroweak observables. Finally, we discuss effects of new particles for muon g-2 on 
the Higgs boson decay /i — )■ 77 and direct productions of these particles at the LHC 
experiment. 



1 Introduction 



The standard model of elementary particles (SM) has been amazingly successful, and cur- 
rently the Large Hadron Collider (LHC) experiment is searching for Higgs boson, which is 
the only particle that has not been observed yet in the framework of the SM. So far, some 
interesting hints of the Higgs boson with a mass of about 125 GeV have been reported at the 
LHC [T] and its discovery may not be in a far future. The Higgs boson with a mass of about 
100 GeV is totally consistent with the electroweak (EW) precision measurements. Therefore, 
the SM well-describes the nature up to the EW scale. 

We, however, think that the SM is not the final theory of elementary particles. The 
hierarchy problem has been a strong driving force to think of the physics beyond the SM, and 
many ideas, such as supersymmetry, extra-dimension, and technicolor, have been proposed 
to solve the problem. The LHC experiment is an ideal place to probe the models for the 
hierarchy problem, and the searches are currently going on. The first round of the LHC 
searches, however, does not show any serious deviations from the SM prediction. It seems 
that this negative search result starts creating a tension between new physics models and 
the requirements for the hierarchy problem. Therefore it will be good time to consider other 
approaches to think about new physics beyond the SM. 

There is a different approach to consider the physics beyond the SM, based on a consider- 
ation for unexplained experimental results. The muon anomalous magnetic moment (muon 
g-2) is one of the most precisely measured observables [2]. The theoretical prediction from 
the SM by several groups has suggested that there is a discrepancy between the experimental 
result and the SM prediction [3j: 

6af, = a;fP - af^ = (26.1 ± 8.0) x 10'^°, (1.1) 

where Sa^ is a discrepancy between the experimental result (o^^p) and the SM prediction 
(a^^). There are many discussions |3lll] on uncertainties of the SM prediction to understand 
whether this anomaly is real or not. If the discrepancy can not be explained by the SM, this 
would be an evidence of the physics beyond the SM. At present, there seems no satisfactory 
explanation for this anomaly within the SM. Therefore, it is worth while considering the 
physics beyond the SM seriously in order to explain the anomaly. 

It is interesting to note that the anomaly of muon g-2 is the same size as 1-loop contribu- 
tions induced by the EW gauge bosons in the SM. This suggests that in order to explain the 
anomaly, new particles with masses of EW scale are required if the interactions are of order 
of EW gauge couplings. 

Although the anomaly of muon g-2 has been discussed in the context of new physics 
models (for example, see Ref. [5] in MSSM, Ref. [6] in extra dimension model, and Ref. [7\ in 
Little Higgs model) motivated mainly by the hierarchy problem, we take a "bottom-up" style 
approach Starting from introducing new interactions with muon to generate new contri- 
bution to muon g-2, we try to capture important features of new physics models, although 

See also Ref. [S^ for the recent related studies. 
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we can only discuss a part of complete models. In this paper, we consider models where 
muon has new Yukawa type interactions with new particles to generate new contributions to 
muon g-2. We find that in a model where the chirarity fiip of muon occurs only in external 
line of muon g-2 diagrams, it is necessary to introduce the relatively large Yukawa coupling 
and the EW scale new particles. On the other hand, in a model where the chirality flip can 
happen in the internal lines of muon g-2 diagrams, we can obtain favorable g-2 contributions 
without large Yukawa couplings, and they arc also consistent with the precision EW observ- 
ables. We also expect these particles may be observed directly and/or indirectly at the LHC 
experiment. 

This paper is organized as follows. In the next section, we consider new physics models 
where muon has new Yukawa couplings to explain the anomaly of muon g-2. One model 
is that only right-handed muon has new Yukawa coupling. Another one is that both right- 
and left- handed muons have new Yukawa interactions to enhance the contributions to the 
muon g-2. We show the consistent region of parameter space. In Sec. 3, we discuss effects 
of the new physics on EW observables and whether the parameter space which is consistent 
with the muon g-2 is also consistent with the precision measurements. In Sec. 4, we discuss 
phenomenology of these new physics models at the LHC. Especially we show the possible 
effect on the Higgs boson decay of h — )■ 77, and the direct production cross sections of new 
particles at the LHC. Sec. 5 is devoted to summary. 

2 New physics models for anomaly of muon g-2 

In order to explain the deviation from the SM prediction of muon g-2 by new physics, muon 
has to have new interactions with some charged particles. Since the effective operator for 
muon g-2 couples to photon, the loop effects via the new interactions of muon with charged 
particles could induce the extra contribution to the muon g-2. 

In this paper, as such new interactions, we consider new Yukawa-type couplings with right- 
handed or/and left-handed muon. We consider two cases: (1) Model where only right-handed 
muon has new Yukawa interaction and (2) Model where both right-handed and left-handed 
muons have new Yukawa interactions. 

2.1 Model where right-handed muon has new Yukawa interaction 
with SU(2) singlet scalar (0) and singlet Dirac fermion (x) 

We consider the following Yukawa interaction where the right-handed muon couples to new 
S\J{2)l singlet fermion x and singlet scalar 0: 

>C = -yNji'RXL(t> - m^XnXL + h.c. - mJ^V H ■ (2-2) 

Here jiR is the right-handed muon. QED charges of new fermion x ^'^^ scalar are and 
= — 1 — Q^, respectively, in order to have the gauge invariant Yukawa interaction. The 
masses of x ^iid are denoted by and m^, respectively. In order to simplify the model. 
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we have assumed a Z2 parity, under which the SM particles are even and new particles 
and X are odd. This type of the parity may be interesting for the dark matter. We also note 
that even if we impose the Z2 parity, the right-handed electron, for example, can have similar 
Yukawa interaction with and x- This will cause severe flavor mixing problem. We will not 
discuss the complete model here, however, we implicitly assume that x or has approximate 
muon flavor number so that the flavor mixing is strongly suppressed. 




Figure 1: Feynman diagrams for muon g — 2 



Since the right-handed muon couples to x and 0, their radiative corrections induce the 
muon g-2, aj^™ as shown in Fig. [H The contribution a"*^"" is given by 



a 



2 2 

new i^jV fi 



167r2 



[QxiCii + C2i)(0, X, -q) - + C22)(0, 0, x; g,p - q)] ,(2-3) 



where p and p — q are momenta of external muons and g is a photon momentum and a 
limit — )■ is taken. Here Cx{A, B,C;pi,p2) {X = 11, 21, 12, 22) are so called Passarino- 
Veltman functions and their definitions in this paper are shown in Appendix. The explicit 
forms are 



, ^ w , s 1 2-3y -6y^ + + 6y\ogy 



(C12 + C22)(0,0,x;9,P- 9) = — 771 u > (2-5) 

where y = m'^/m'^, = and the higher order terms of 0{rn?^/rn?^) are neglected. The first 
term in Eq. (12. 3p comes from Fig. [TJ^a) and the second is from Fig. [H^b). 

In Fig. El we show the contribution to a^"" as a function of Q^^ and m^. Here we take y^ = 
2.5 and = 200 GeV. We show contours for (a;i'="/10-^°) = 2.1, 10.1, 18.1, 26.1, 34.1, 42.1 
and 50.1, from right to left, corresponding to —3(7, —2a, —lex, Oct, la, 2cr and 3(7 deviation 
from the measured value, respectively. We note that y^v dependence of a°*^^ is trivial, that 
is, a^'^^ is proportional to yj^, as shown in Eq. ( 12. 3p . For example, if one takes yN = 1, the 
values of a"^"" in Fig. [2] reduce by a factor (1/2.5)^ = 0.16. 
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Figure 2: New physics contribution to muon g-2 (a^^™) as a function of and m^. 
Here we take = 2.5 and = 200 GeV. We show contours for {a^^^^ /lO^^^) = 
2.1, 10.1, 18.1, 26.1, 34.1, 42.1 and 50.1, from right to left, corresponding to —3a, — 
2cr, — Icr, Oo", la, 2a and 3a deviation from the measured value, respectively. 



One can see that the region with > — 1 (which corresponds to < 0) is disfavored 
by the data of muon g-2. It is interesting to notice that the neutral fermion {Q^ = 0, 
corresponding to = —1) is difficult to explain the anomaly in any values of m^, and 
on the other hand, the neutral scalar {Q^ = 0, corresponding to = —1) can potentially 
accommodate the anomaly if the scalar is not so heavy. It is also interesting to note 
that the multi-charged fermion and scalar (such as = —2, — 3, ■ ■ ■ , corresponding to 
Qff, = 1, 2, ■ ■ ■ ) are also favored by the muon g-2 anomaly if is a right value for the 
anomaly. Therefore, the anomaly of muon g-2 constrains the QED charges as well as the 
mass scale of new particles. 

In Fig.[3l we show a^^^ as a function of and for (a) Q^^ = — 1 (corresponding to = 
0) and un = 2.5 and (b) = —5 (corresponding to Q^p = 4) and = 1.5. Contours for 
(a^'^^/lO-^") = 2.1, 10.1, 18.1, 26.1, 34.1, 42.1 and 50.1 (corresponding to -3a, - 2a, - 



la, Oa, la, 2a and 3a deviation from the measured value) are shown, similar to Fig. 121^^1 As 



one can see, new particles with their masses of about 100 GeV and the coupling ~ 0(1) 



If one changes the value of Yukawa couphng jjn, one should simply multiply the value of aji°* in Fig. |3] 
(a) and (b) by a factor (y|)^ and (f^)^, respectively. 
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Figure 3: New physics contribution to muon g— 2 [a^^^) as a function of and for 
(a) = —1 and and yj\f = 2.5 (b) = —5 and = 1.5, respectively. Contours for 
(a;i^™/10-^°) = 2.1, 10.1, 18.1, 26.1, 34.1, 42.1 and 50.1 (corresponding to -3a, -2a, - 
la, Oct, la, 2a and 3a deviation from the measured value) are shown, similar to Fig. |5J 



are strongly favored by the anomaly of muon g-2 when = — 1- When the QED charge of 
X decreases further negatively (from —1 to —5), the muon g-2 gets larger even if the Yukawa 
coupling gets smaller, as shown in Fig. [3] (b). 

Since the anomaly of muon g-2 requires relatively light new particles (~ 0(100 GeV)) 
and/or relatively large Yukawa coupling (~ 0(1)) in this scenario, the EW precision observ- 
ables can be affected by these new particles and new interactions. Therefore, we will check 
the effects of these particles on the EW observables in a next section. 

We would also hke to point out that the anomaly of muon g-2 suggests the existence of 
relatively light new particles and/or multi-charged particles. Therefore, it is very interesting 
to know whether these particles can be found directly or indirectly at the LHC We will study 
the direct productions of these particles and the effects on Higgs decay to 77 at the LHC in 
a later section. 

2.2 Model where both right- and left-handed muons have new 
Yukawa couplings 

Unlike the previous case, here we consider that both right- and left-handed muons have new 
Yukawa interactions: 

C = -yLL2^XR-yR^R<PXL-m^XLXR + ^-c., (2.6) 
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where L2 (= {i^^il, f^L)^), ^ (= (01; 02)'^) and are the second generation SU(2) doublet 
lepton, SU(2) doublet and singlet scalars, respectively, and x is SU(2) singlet fermion, whose 
mass is m^. Here we have assumed the Z2 parity again in order to simplify the model. 
Under the Z2 parity, the SM particles are even, and new particles 0, $ and x are odd. We 
also implicitly assume the approximate muon flavor symmetry, under which x or (0, $) has 
muon flavor number, so that the flavor mixing of this type of new Yukawa couplings are 
strongly suppressed. The model we discuss in this paper may be a part of the cornplete 
model. However, we think this part of the complete model is crucial for the muon g-2j^i 

The QED charges of new particles are represented in term of QED charge of x field (Qx) 
as follows: 



Q(0i) 


= Qi = 




Q(02) 


= Q2 = 




Q<t> 


= -1 - 


Qx = 



(2.7) 
(2.8) 
(2.9) 

Since = Q2, (f) and 02 can mix each other. For example, the following gauge invariant 
term induces the 02 — mixing mass terms after the electroweak symmetry is spontaneously 
broken: 



C = -AM(i7"f$0"f) + h.c. 



XMv 



+ 



(2.10) 



Here we parameterize the mass terms for and 02 as follows: 



#4 



£ = - 



m 



11 



m 



12 



''12 "^22 J \ V2 

Thus diagonalizing this mass matrix, we define the mass eigenstates Sj {i 



(2.11) 



1,2) as 



V--S- 



(2.12) 



Here Vij is a mixing unitary matrix, which diagonalizes the mass matrix shown in Eq. (12.111) . 
The mass eigen values are taken to be m^^ < rris^- This type of mixing terms are important 



summary, we briefly comment on the possible completions of this type of model. 
"^^When Qif, = 0, mass term such as — m^002 is also possible. Here we neglect such a mass term for 
simplicity. Even if we include such a mass term, our result does not change qualitatively. 
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to induce large contribution to muon g-2. The contributions to muon g-2 are summarized by 



= -ll/ E { l^2.r + yl\Vu\') (Cn + X, x; P, -q) 

i 

Tfl 



-rri I 

+2?/i?/R^Re(\/2il^i*)Ci2(s„s,,x;g,p-g) , (2.13) 

J 



where p and p — q are momenta of external muons, and g is a photon momentum and a limit 
— )■ is taken. Here we show the explicit expressions of above Passarino-Veltman functions: 



(Cii + C2i)(si,x,X;P,-g) = ^ Q{l-y^^ ' ^ ^ ^ 



c^n(..x,x;P,-g) = , (2.15) 

(C22 + C,2)(....x;.,P-.) = J_l^%±|[±^|z%[i!^, (2.16) 

mi^ 6(1 - ?/i)4 

1 ^-y'i + 2y » In 

"^l 2(1 - 



<^i2(si,Si,x;g,p- g) = — —^ — r\i^ — ' (2-1') 



where i/j = m^/ml^ and we neglected the higher order of 0{rn?^/m'l^. 
The effective operator which expresses the muon g-2 is written by 

C = ^liR^y^" IIlF^u + h.c. , (2.18) 

where f is a vacuum expectation value of Higgs boson, and F^i, is a field strength of photon 
field, and A is a typical scale related to new physics. As one see, the chirality of muon has 
to flip in this interaction. In the case where only right-handed muon has the new Yukawa 
interaction, the chirality flipping of muon happens in the external muon line in the loop 
diagram. On the other hands, in the present case where both right- and left-handed muons 
have the new Yukawa interactions, the chirality flipping can occur in the internal fermion 
line, which is proportional to the mass of the fermion x- That is why there are terms which 
are proportional to yRyLfn^/m^ in Eq. f l2.13p . This contribution enhances the effects of muon 
g-2. Therefore, it is very important to explain the anomaly of muon g-2. 

In Fig. m we show a^"^ as a function of and (which is a mass of lighter scalar 
state) in case of (a) = -1 and (b) = 0. Contours for (a°^"/10-i°) = 2.1, 10.1, 
18.1, 26.1, 34.1, 42.1 and 50.1 (corresponding to —3a, — 2a, — la, Oa, la, 2a and 3a 
deviation from the measured value) are shown. Here we took mf^ = 171^2, ^^12 = (50 GeV)^, 
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Figure 4: New physics contribution to muon g-2 {a^'^'") as a function of and rris^ in case of 
(a) = -1 and (b) = 0. Contours for (a;i^"/10-i°) = 2.1, 10.1, 18.1, 26.1, 34.1, 42.1 
and 50.1 (corresponding to — 3o", — 2a, — la, Oct, la, 2a and 3a deviation from the 
measured value) are shown. Here we assume that ml-^ = m\2i ^^12 = (50 GeV)^. We also 
take (a) yi = -ya = 0.4 and (b) yL = yR = 0.4. 

and (a) yL = —yn = 0.4 (b) yi = Vr = 0.4. Note that in order to generate the positive 
contribution to muon g-2 in Fig. H] (a), the sign of new Yukawa couplings had to be taken as 
VlVr < 0. On the other hand, in Fig. H] (b), the sign of yiVR should be positive. As shown 
in Fig. m compared to Fig. [3], even smaller Yukawa couplings and heavier new particles can 
accommodate the anomaly of muon g-2 because of the enhancement mentioned above. As 
can be seen from Fig. HJ when |/l ~ ~ 0(1), the EW scale new particles are expected in 
order to explain the anomaly of muon g-2. Therefore, we will analyze the effects on the EW 
precision observables and show the consistent region of parameters in a next section. 

3 Effects on electroweak observables 

In the previous section, we showed that the relatively light new particles are required in order 
to explain the anomaly of muon g-2. In addition, in the case where the only right-handed 
muon has new Yukawa interaction, the new Yukawa coupling should be relatively large for 
the anomaly of muon g-2. Therefore, we should check if the relatively light particles and 
the relatively large Yukawa coupling are consistent with the EW precision measurements. 
Although models which we consider in this paper may be only a part of complete models, we 
should know how these new particles affect EW observables. 

In this section, we adopt the formalism in Refs. [TOl [TTl [12] in order to include the oblique 
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corrections as well as the vertex corrections in the EW observables. First we briefly summarize 
the formalism pB[m[T2]. 

In the presence of EW to TeV scale physics, it is well known that the oblique corrections 
in gauge boson self-energy are important. They are parameterized by Peskin-Takeuchi's S, 
T, U parameters [T3| [HI [T5] : 



aS Uzz{Ml)-Uzz{0) c^w ^z^{Ml) Ii,,{M\ 



2- 



isl^cl Ml cwsw Ml Mf 



(3.19) 



z 



Uww{0) UzziO) 



= (3.20) 



Ml 

aU ^ UwwjM^) - nw(0) 2 ^zzjMl) - UzzjO) 
^s^w Ml 

Uz,{Ml) n,,(M|) 

-^swCw^H ^2 • (3-21) 

Here we use the notation cw and sw to refer to the cosine and sine of the weak mixing angle 
and C2W = ~ In addition, in Refs. [TOt [TT] . Rz and Rw parameters are introduced in 
order to account for the smaller corrections: 



(xRz dUzzip'^ 



HzziMl) - UzziO) 



p^=Ml Ml 



aRw IlwwiMl) - Uww{M^) ^ww{M^) - ^ww{0) 



islr Ml - M^ M^ 



(3.22) 
(3.23) 



Furthermore, new particles have effects on the running QED coupling constant q;(M|) [31IT2]: 

a(M|) = J- (3.24) 

1 - Aaiep(M|) - Aa;fj,(M|) - Aatop(M|) - A«„ew(M|) 

Here Aaiep(M|), Aa'^J^{Ml), Aatop{Ml), and Aanew(^l) are the leptonic, the five-flavor 
hadronic, the top quark and the new physics contributions to the running of the QED coupling 
constant respectively. The new physics contributions to the running of the QED coupling 
constant Aanew{Ml) are defined as 



Aa„ew(M| 



H,,(M|) U,,{p'^ 



Ml p2 



(3.25) 

p2=0 



These oblique parameters contribute to EW observables, as shown in Refs. [TOl [TTl [T2] . 



If muon has 0(1) new Yukawa couplings]^ it induces non-universal vertex corrections to 



the Zfj,~^fj, coupling. The standard model coupling of Zfi^fj, is given by 



^^7,[9r''PL + 9T''PB], (3.26) 

Cw 



'^^In this paper, we assume that electron and tau do not have 0(1) extra Yukawa couphngs, for simphcity. 
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Figure 5: (a) of the EW observables listed in Table. [T] as a function of the Higgs boson 
mass m/j within the SM. (b) contours shown in dashed lines in S-T plane, assuming that 
other oblique corrections (f/, Rw) Rz) and vertex corrections are zero. Reference Higgs 
boson mass and top quark mass are taken to be 125 GeV and 173.2 GeV, respectively. We 
also show the predicted S-T values (shown in solid lines) in the SM, varying the Higgs boson 
mass and top quark mass from the reference values. 



where g^^^^ are the standard model couplings, and their tree level contributions are g^^'^ = 

— I + and g^''^ = s^. When we take into account of corrections to Z vertex via 
the new Yukawa coupling (including the corresponding wave function renormalization) , the 
Z/i+yU" coupling is modified by 



. g 

cw 



Here the vertex corrections generated by new particles are parameterized by ^g^ j^. 

The corrections to W iiv^ vertex generates the corrections to the Fermi constant Gp and 
therefore, we parameterize it as A^g*. The A^g for the new physics contributions from the 
vertex and box diagrams to the /x-decay process is defined as 

Here Gf^+°^- is the muon decay constant including effects of SM radiative corrections and 
new physics oblique corrections except the vertex and box corrections from new physics. 

Using the formalism in Ref. [T0| [TT] . one can calculate EW observables listed in Tabled] 
without assuming lepton universality. First, we show a result of the SM fit of the EW 
observables in Fig. E] (a) as a function of Higgs boson mass rrih- The best fit point is at 
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Table 1: Experimental data and theoretical predictions of electroweak observables. The 
predicted values of the SM and sample model discussed in Sec. 3.2 are shown. 
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rrih = 91 GeV. This is consistent with the result shown by the LEP electroweak working 
group [T8|. In Table [H we show the SM fit when Higgs boson mass is assumed to be 125 
GeV, which has been suggested by the latest LHC data. In Fig. [5](b), contours are shown 
in dashed lines in S-T plane. Here we assume that other oblique corrections {U, Rw, Rz) 
and vertex corrections are zero. Reference Higgs boson mass and top quark mass are taken 
to be 125 GeV, 173.2 GeV, respectively. The predicted values of S-T parameters are also 
shown by varying the Higgs boson mass and top quark mass from the reference values. We 
note that the SM fit with lighter Higgs boson is good except for the muon g-2, and small 
positive S and T (S" ~ 0.05, T ~ 0.1) can decrease the further. Therefore, we need new 
physics which largely contributes to muon g-2, but whose effects on other EW observables 
are small. In this section, we analyze the EW observables in models discussed in the previous 
section. In this analysis, we assume irih = 125 GeV. 



3.1 Model with SU(2) singlet scalar (0) and singlet Dirac fermion 

(x) 

Since the new scalar (f) and new fermion x are SU(2) singlet, they do not couple to W boson. 
But if it has a QED charge, it can couple to photon and Z boson. Here we list the corrections 
to self-energy functions of gauge bosons in this model. 
The singlet fermion (x) contributions are given by 

ni^l^(/) = 

4"1(/) = -^^K^o(x,x)-/{i?i(x,x) + i?2i(x,x)}-2(l-e)i?22(x,x)], 
n^^Hp') = -^Ql[mlBo{xx)-p'{B^{xx) + B2lix,x)}-2{l-e)B22{x,x)], 

nW(p2) ^ 9_^^^^2^Boixx)~p'{Biix,x) + B2i{xx)}-2{l-e)B,,ixx)]- 

(3.29) 

Here Bx{i,j) = Bximf^rn^-^p) {X = 0, 1, 21, 22), which are Passarino-Veltman functions 
and are shown explicitly in Appendix. We use the dimensional regularization in space-time 
dimension D = 4 — 2e. 

The singlet scalar (0) contributions are given by 

niJW(p') = 

^fzip') = -^Qlstv[B22{<l>,<l>)-lA{<p) 

n?,H/) = ^QI[b22{<PA)-Ia{<p)Y 

^^,%') = -4^Q|^'v.{i?22(0,0)-^A(0)|. (3.30) 
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Here Passarino-Veltman functions i?22(^,j)(= B22{Tnf,m'j;p)) and A{(f)) are shown in Ap- 
pendix. We can easily show that Peskin-Takeuchi STU parameters in this case vanish 
{S = T = U = 0), since the new particles do not have SU(2)l interactions. Thus, the 
leading contributions to the oblique corrections are Rz parameter and Aancw(^l): 




(3.31) 



.(Ml) 



9tt 



97r 



1 + 



12ml 
5M| 





1 tan 



1 tan 



-1 



4mj 



4m2 _ 
Ml" \ 



(3.32) 



where we assume 2m^ > and 2m0 > Mz- 

As discussed in the previous section, in order to explain the anomaly of muon g-2, the 
new Yukawa coupling in this model should be relatively large. In this case, the potentially 
large vertex corrections may be expected, as shown in Figj6l The results are expressed by 
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(3.33) 



Vn 
167r2 



+Qx^w 1^ - 2C24 - M|(Ci2 + C23) + m^Coj (x, 0, x; g - p,p) 

-s^(i?o + 5i)(0,X;p)], (3.34) 



where p, q—p and g are muon, anti-muon and Z-boson momenta, respectively. Here Passarino- 
Veltman functions Cx {X = 0, 12, 23, 24) and Bx {X = 0, 1) are explicitly shown in 
Appendix. 
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Figure 7: (a) Effect of vertex corrections on i?^. = — R'^ are shown as a function 
of and in case of = —1 and yx = 2.5, where R^ includes all corrections, and R'^ 
contains all corrections except for the vertex corrections, (b) as a function of and 
in case of = —1 and yx = 2.5. 

The vertex corrections affect EW observables such as i?^, Ap^, and A^. Among these 
observables, we find that the effect on the R^ is the most important in the fit because i?^ is 
most precisely measured. In FigJ7| (a), we show the effect of vertex corrections in i?^, that 
is, SRfi = Rfj, — R'^ as a. function of and in case of = —1 and yx = 2.5. Here 
Rfj^ is a theoretical prediction in this model containing all corrections, on the other hand, 
R'^ contains all corrections except for the vertex corrections. The difference shows the 
effect of the vertex corrections. Note that the size of experimental error of i?^ is 0.033 for 
la, as shown in Table. [TJ As one can see in Fig. El^a), the vertex correction can change the 
prediction of i?^ about la in the region where the muon g-2 can be explained by the new 
physics contributions. Since the SM prediction of i?^ is smaller than that of the experimental 



14 



500 



400 



-5 vn = 1-5 



^ 300 % 



200 r 



100 



' I I I IV ' ' ' I ' ' ' ' I ' 'V ' I ' ' ' ' I ' ' I ' ' ' ' I ' ' ' I. 



1-0.01 



-0.02 



\ -0.05 



-0.006 / 
/ 
/ 
/ 

/ - 
/ 

/ 

/ 

/ 

/ 

/ 



/ /-0.004 



> 



500 



400 



300 



200 



Qx = ~5 ?/Ar = 1.5 



100 



40 



I I |1 I 1 I I 1 I 1 I 1 J I I I I I I I I I I I I I I I 



35 /30 
//29 



''/ - 

mf. ; 

'■\%\\ 125 



100 



200 



300 



400 



500 



100 



(a) 



[GeV] 



200 300 400 500 

[GeV] 

(b) 



i?' is shown function of 



Figure 8: (a) Effect of vertex corrections on R^. bR^ = R^ 

and in case of Qy. = —5 and i/n = 1-5, similar to Fig. [Tl^a). (b) a function of 

and in case of = —5 and yjy = 1.5. 



result (shown in Table. [T]), the vertex corrections do not help to improve the prediction of 
R^. Therefore, the region of small x mass is very constrained. 

In Fig. [7] (b), is shown as a function of and in case of = — 1 and z/n = 2.5. 
As one can see from the figure, is large in the region of small because of the large vertex 
corrections, discussed above. In the region of larger and m^, it is difficult to explain the 
anomaly of muon g-2 because new particles are too heavy, and hence the gets larger. The 
minimum of x^ is around ~ 200 GeV and ~ 100 GeV, and we see that the light new 
particles are favored by the EW observables including muon g-2. 

In Fig. M (a), the effect of vertex corrections on R^, 6R^ is shown in case of = —5 
and yN = 1-5, similar to Fig. [7] (a). The vertex corrections can be as large as Icr of the 
experimental error of i?^, and they increase the x^- In Fig- El (b), the x^ is shown as a 
function of and m^. As one can see, the region of small rUy. is disfavored by the vertex 
corrections of i?^, as shown in Fig. [8] (a). The minimum of the x^ is around rriy. ~ 300 GeV 
and ~ 100 GeV. Therefore, the EW precision measurements (including muon g-2) prefer 
the relatively light new particles in this case. 
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3.2 Model with SU(2) doublet ($) and singlet (0) scalar, and sin- 
glet fermion (%) 

The singlet fermion contributions are same as those in the previous case. Here, we hst the 
one loop contributions to gauge boson self-energies induced by scalar sectors. 



3_ 



(3.35) 



47r 



(3.36) 



ng(p^) 



47r2 



g?522(0i>i)- 7.^(01) 



■i ^ 



_e 9_ 

47r^ 



Ql ( ^ - Qisl\ |i^22(01>l) - ^A(0i) 



X \ B22{Si,Si) - -A{Si 



(3.37) 



(3.38) 



where Passarino-Veltman functions Bx{i,j){= Bximf^mf^p)) and A{i) are shown in Ap- 
pendix. Unhke the previous model, SU(2) scalar doublet can contribute to the STU param- 
eters. Therefore, the dominant quantum corrections are represented by S and T parameters. 
To understand a behavior of S and T parameters, here we show the approximate expression 
of S and T parameters. 



S ^ + 

.2 



T ~ 



(A)^ + 



(3.39) 



16 



Here, for simplicity, we assume that V22 = Vu = 1, Vu = V21 = so that S2 state originates 
from 02 component of SU(2) scalar doublet $, and A = {m'^^ — fn'^J/m'^_^ where m^^ = "^82 
in this case. Here we ignore the higher order terms of A and m\/m?^^ in Eqs. f l3.39p . Note 
that A parameterizes non-degeneracy in the SU(2) scalar doublet is a hypercharge of 

$, Y$ = \ + Q2- We note that the non-degeneracy of SU(2) scalar doublet generates non-zero 
T and non-zero S, and T is always positive, but S can be positive or negative depending on 
signs of Y$ and A in this model. 
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Figure 9: contours (without including muon g-2 result) in S-T plane ( = 22 to 34 
in dashed lines) and S-T values (solid lines with points) in this model. Here we take ?/l = 
—Vr = 0.4, = —1 (corresponding to Qi = 1 and Q2 = 0), and mj^ = = (300 GeV)^. 
Two solid lines correspond to two different values of T/ix2) ^12 ~ 

(50 GeVf and (250 GeV)^, 

varying m^^ from 200 GeV to 400 GeV at 50 GeV step. 

In Fig. [HI we show numerical result of contours in S-T plane (without including muon 
g-2 result), shown in dashed lines. To draw the contours in S-T plane, we assume that 
other oblique corrections {U, Rz, and Rw) and vertex corrections are zero. As one can 
see from the figure, slightly positive S and T (5* ~ 0.05 and T ~ 0.1) are favored by EW 
observables. In Fig. [HI we also show the predicted values of S-T parameters in this model, 
shown in solid lines. Here we took = —yn = 0.4, = —1 (corresponding to Qi = 1 and 
Q2 = 0), and m^^ = = (300 GeV)^. Two solid lines correspond to two different values of 
m^2) ^12 = (50 GeV)^ and (250 GeV)^. The points on the sohd lines represent the predicted 
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values of S and T, varying m^^ from 200 GeV to 400 GeV at 50 GeV step. 



When mfi 



m: 
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(300 GeV)2 and mj^ = (50 GeV)^, m,, 



296 GeV and ~ 



304 

GeV. As roughly shown in Eqs. (I3.39p . when m^^ = 200 GeV and (^2 = 0, A > and 
y$ = ^ > 0, and hence S > and T > 0. As rriip-^ becomes larger and about 300 GeV, both 
5* and T get closer to zero because the doublet scalars are almost degenerate. Then when 
m<^j is increased further, S gets negative because A < 0, but T > 0. In Fig. [9l one can see 
this behavior. Even when we increase mf2, we can see that this behavior is almost same. 
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Figure 10: Same as Fig. [9] except for Ul = Vr = 0.4 and = (corresponding to Qi = 
and Q2 = -1). 

In Fig. [TOl the same figure is shown except for i/l = = OA and = (corresponding 
to Qi = and Q2 = -1). When Q^ = {Q2 = -1) and m^^ = 200 GeV, F$ = -i < and 
A > 0, so that S is negative. As m^-^ gets larger, S gets larger. For ~ 300 GeV, S' ~ 
because the doublet scalars are almost degenerate. For m,^^ > 300 GeV, S becomes positive. 
This behavior is different from the previous case with Q2 = since the sign of is different. 
Note that T is always positive unless SU(2) scalar doublet is degenerate. As one can see in 
both cases, the small non- degeneracy of SU(2) scalar doublet can make better. 

In our numerical analysis, we also include the vertex corrections in Zn~^fi~, Zv^j^v^ and 
WfiUfj, vertices even though the Yukawa couphngs i/l and yn can be small in this model. The 
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expression of the vertex corrections is summarized as follows: 



167r2 



+ E ^x4l^2.r 1^ - 2^24 - M|(Ci2 + C23) + mJCoj (x, s„x;q- P,p) 



(3.40) 



Q<^Sw 1 ^2*% - Qr^sl^V^^ij \ ViiV*jC24{Si, X, Sj]P, Q - p) 



+ Yl QxSw\Vu\' 1^ - 2C24 - MliCu + C23) + mJCoj (x, s^, X; 9 - 



^9T 



16n 



J2sw\Vii\''iBo + B,){si,x;p) 

C24(0i,x,0i;9 -p,p) 



(3.41) 



1 



2 ( 2 ^wQfki 



+QxSw 1^ - 2C^24 - M|(Ci2 + C23) + m^Coj (x, 4>i, x;q- P,P) 



--(fio + 5i)(0i,x;p) 



(3.42) 



where p and g are momenta for muon in Ag'^ ^(muon neutrino in Ag'^^) and Z-boson, respec- 
tively. Correction to /zz/^M^-vertex is written by 



G 



87r2 



^ |^2ipC'24(0i,Si,x; -g,p) 



i J (5o + 5i)(0i, x; P - g) + 5^ I V2i|'(5o + B^){si, x; p) 



(3.43) 



where p, p — q and g are momenta for muon, muon-neutrino, and W-boson, respectively. Here 
we use approximation = 0. 

In Fig. [m we show x^ contours as a function of and rris^ in case of (a) = —1 and 
Ul = —Ur = 0.4 and (b) = and i/l = Vr = 0.4. Here we have taken = (50 GeV)^, 
and mil = ''^22 ■ ^<t>i chosen so that the x^ becomes the smallest, provided that m^-^ > 
105 GeV for = —1 and m^^ > 46 GeV for = respectively. In both cases shown in 
Fig. [TT], the x^ can be as small as 22, which is smaller than those in the previous case where 
only right-handed muon has new Yukawa coupling. As seen from Fig. [HI the relatively light 
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Figure 11: contours as a function of and m^^ in case of (a) = — 1 and 
0.4 and (b) = and i/l = Vr = 0.4. Here we take mf-^ = and mfg = (50 GeV)^. 



X and Si are favored by the EW observables. In Table. 1, we show the predicted values of 
the EW observables in case of = 0, m^.^ = 300 GeV, rriy^ = 200 GeV, m^^ = ml^ = (250)^ 
GeV^, and 

772,^2 — (50)"^ GeV^. As shown in Table. 1, the fit is quite good. Therefore, it will 
be important to know what kind of effects we expect from these particles at the LHC. In a 
next section, we will discuss the possible effects of these new particles at the LHC. 

4 Phenomenology at the LHC 

In the previous section, we show that in order to explain the anomaly of muon g-2, it is 
strongly suggested that there should be the EW scale new particles. It is also suggested 
that multi-charged particles may be favored by the anomaly of muon g-2. Therefore, these 
particles may be reachable directly and/or indirectly at the LHC. Here we discuss the effects 
of these new particles in the Higgs decay /;.—)■ 77 and their direct productions at the LHC. 

4.1 /i — 77 

The SU(2) singlet and doublet scalars can couple to the Higgs boson through the following 
scalar interactions: 

-K4M { {H^^)(j)'^ +h.c. } , (4.44) 

where if is a standard model Higgs boson doublet. In case where only SU(2) singlet scalar 
(j) exists, only coupling k,i is non-zero, but in case where both SU(2) singlet (j) and doublet $ 
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scalars are there, all couplings are alloweccj. 

In the case where only right-handed muon has new Yukawa interaction, there is SU(2) 
singlet 0, and its interaction with Higgs boson is 

£ = -Xvh(f)^(f), (4.45) 

where A = ki in this model and f is a vacuum expectation value of Higgs field, v ~ 246 GeV. 
When Q^, which is the QED charge of 0, is non-zero, can contribute to the Higgs decay 
h — )• 77. The decay width r(/i — )■ 77) is given by 

2 s 

^^^^^^^ = ^^1^^"^^)!'' (^-^^^ 

where the amplitude S{mh) at leading order is written by 



,2 



5(m,) = -Ft{n) - Fwirw) + Q%\^F^{t^). (4.47) 
Here the first, the second, and the third terms are top, W, and contributions, respectively, 

2 

and Tx = {x = t, W, 0). The function Ft, Fw, F^ are given by 

Ftir) = r-i[l + (l-r-i)/(r)], (4.48) 
Fw{t) = 2 + 3r-i + 3r-i(2-r-i)/(r), (4.49) 
F^{t^) = r-i[-l + r-V(r)], (4.50) 

where /(r) = arcsin^(y^) for r < 1. Branching ratio of /i — )■ 77 is approximately given by 



BR{h ^ 77) T{h -> 77) 



BR{h 77)sM T{h 77) 



SM 



\Ft{n) - FwM + QlX^F^ir, 

4> 



\Ft{rt) - Fw{tw) 



, (4.51) 



because the total decay width does not change much, even if r{h — 77) changes. Here 
BR{h — !■ 77)sM and T{h 77)sm are the SM predictions at the leading order. Notice that 
this ratio is a function of and Q^A when Higgs boson mass mh is fixed. As we discussed 
earlier, in order to explain the anomaly of muon g-2, larger may be favored. As can be 
seen in Eq. fl4.5ip . larger induces larger effect on /i — )■ 77. 

In Fig. [121 the ratio T{h — i- ■yj)/T{h — )• 77)sm is shown as a function of and Q^X- Here 
we took nih = 125 GeV. As one can see from the figure, the effect can be significant if Q^A 
is large and m<^ is O(IOO) GeV. For example, as we discussed, when = 4, A = 1 (— 1) and 
rrifp = 200 GeV, T{h — )■ 'yy)/T{h — )■ 77)sm = 0.12 (2.7). Unfortunately, we can not predict 



When the QED charge of cj) {Q4,) is zero, terms such as {H'^^)(j) and are also possible. In our 

analysis, we assume that such terms are small, for simplicity. Even if we include such terms, our result does 
not change qualitatively. 
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Figure 12: Higgs decay width /i — )■ 77 normalized by the SM prediction f^j^^^^-^^_^ 
function of and Q'iX. Here we take ruh = 125 GeV. 



as a 



the branching ratio BR(/i — )■ 77) in this model because A is unknown. However, it is not 
surprising that BR(/i — )■ 77) is very different from the SM prediction. 

In the case where both right- and left-handed muons have new Yukawa interactions, the 
SU(2) singlet and doublet scalars have interactions with Higgs boson, as shown in Eq. f l4.44p 
and they generates the following interactions: 



£ = — K2vh(j)[(j)i — Xijvhs*Sj, 



(4.52) 



where 



M 



They contribute to the decay width of /i — )■ 77, 



(4.53) 



BR{h 77) 
BR{h -> 77)sM 
\FtiTt) - Fwirw) + QIk2^F^{t^,) + Ql A 



(4.54) 



When = —1, which corresponds to Qi = 1 and Q2 = (as an example we discussed 
above), only 0i can contribute to h 77, because Si does not couple to photon. The 
result is same as one in Fig. [T2] for X = k,2 and Q^p = Qi = 1. In this case, the effect 
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can be significant if K2 is large. For example, when m^^ = 250 GeV and K2 = 1 and — 1, 
T{h — )• 77)/r(/i — )• 77)sM = 0.95 and 1.05, respectively. 

When = 0, which corresponds to Qi = and Q2 = —1, scalars Si {i = 1,2) can 
contribute to h 77. Similarly to the previous case, if the Higgs couplings An and A22 are 
large, the effect on h can be significant. For example, when mli = and M = 

we obtain An = n and A22 = 2/? for = k > (i = 1 — 4), on the other hand. An = 2k, 
and A22 = 1^ for = k < (i = 1 — 4). As a result, for m^-^ = m22 = (300 GeV)^, 
T{h -> 77)/r(/i -> 77)sM = 0.90 and 1.14 for k = 1 and -1 (0.71 and 1.56 for k = 2 and -2 
), respectively. Although it is difficult to predict the decay branching ratio of — )• 77 in this 
model because it depends on many unknown parameters in scalar interactions, it is possible 
to have large effect on /i — > 77 process. 



4.2 Direct productions at the LHC 

As we have discussed in the previous sections, the anomaly of muon g-2 strongly suggests 
the existence of EW scale new particles. Models we discussed in this paper may not be 
complete, namely, there may be other new particles in the complete models. Therefore, in 
this paper, we can not discuss the well-defined signatures at the LHC. However, we think 
that it is important to know how many events of new particles can be generated at the LHC. 
Here we present the production cross sections of new particles at the LHC. 
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Figure 13: Pair production cross sections of (a) x case of = — 1: pp {pp) Z/'-f ^ 
X'^X": and (b) 02 in case of Q2 = —1: pp {pp) — )■ Z/j 4>24>2 ^-t the LHC (Tevatron), as 
a function of and m^^^ respectively. Here we take a/s = 8 TeV and 14 TeV for LHC and 
1.96 TeV for Tevatron. 

In the case where the right-handed muon has new Yukawa coupling, singlet scalar (0) and 



23 



singlet fermion (x) exist in the model. When = — 1, only x couples to the standard model 
gauge bosons {Z and 7) because = 0. Therefore, the fermion x can be pair-produced 
via EW gauge interactions. We calculate the pair production cross section of x at the LHC 
{y/s = 8 TeV and 14 TeV) and at Tevatron (a/s = 1.96 TeV) as a function of m^, as shown 
in Fig. [TSTa). To calculate the production cross sections, we have used the MadGraph [19]. 
In order to explain the anomaly of muon g-2, the mass of the fermion x should be of order 
of about 100 GeV. From Fig. [TST a). the production cross section of x is in the range of 
0.5-0.0005 pb (1-0.001 pb) for = 100 - 500 GeV and = 8 TeV (vs = 14 TeV) at the 
LHC and also 0.1-0.001 pb for = 100 - 250 GeV and ^ = 1.96 TeV at the Tevatron. 
If we consider a luminosity of 10 fb~^ for y/s = 8 TeV and 14 TeV, the number of events 
produced at the LHC is of order of 5-5000 and 10-10000, respectively. Therefore, number of 
signal events can be significant. In this case, if is lighter than x can decay to fi and 0. 
If is a stable neutral particle, the final signature may be a fi'^fj,^+ missing energy. This 
will suffer from the standard model background such as W~^W~{-^ fJ'^f^ ) production whose 
production rate is about 0.29 (0.51) pb for y/s = 8 TeV (14 TeV) at the LHC. Therefore, the 
detail study will be very important to discuss the discovery potential of such a signal event. 
However, the detail signature will depend on the particle spectrum in the complete model. 
Thus we will leave the detail study. If we consider multi-charged the cross section increases 
by a factor Q^, and the cross section further becomes larger, and hence the constraints from 
the LHC will be very important. 1^ 

In the second model where one has a SU(2)l doublet scalar in addition to and the 
scalars can be pair-produced via EW gauge interaction. For example. Fig. [13] (b) shows the 
production cross sections for the charged 0^02 scalars as a function of m^^ in case of = 0, 
which corresponds to Qi = and Q2 = —1. We show the predicted values for y/s = 8 TeV 
and 14 TeV at the LHC and ^/s = 1.96 TeV at the Tevatron. As can be seen from the figure, 
the cross sections are in the range of 0.1-0.0001 pb (0.2-0.0005 pb) for m^^ = 100 — 500 GeV 
at = 8 TeV (^i = 14 TeV) and also 0.02-0.0001 pb for for m^^ = 100 - 250 GeV at the 
Tevatron. For simplicity, we have assumed that there are no 0-02 scalar mixing. In this case, 
01 -f- 02 are also produced via ly-boson. If masses of SU{2) ^ doublet scalar are degenerated, 
we have checked that the cross section is similar to one in Fig. [12] (b). In general, the signal 
events have the form pp {pp) V ^ 'Pl'Pj where V = Z^'y^W and z,j = 1,2. Therefore, 
Fig. [13] shows the typical values of the cross section in this model. 

Although we can not discuss the complete signatures because our sample models are not 
complete, we stress that the searches for the production processes of these new particles 
via EW interactions are very important in models where the anomaly of muon g-2 can be 
explained by these new particles. 

'^''in addition, pair production is also not negligible in this case. 
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5 Summary 



The current LHC experiment does not show any deviation from the SM predictions. This 
result puts severe constraints on new physics models. Especially, colored particles are strongly 
constrained since the LHC as a hadron collider can easily produce the colored particles 
via QCD interactions. For example, in the constrained minimal supersymmetric standard 
model (CMSSM), gluino as well as the first and the second generation squarks are strongly 
constrained and their typical lower mass limits are about 1 TeV [20]. Therefore, models 
motivated by the hierarchy problem are getting more constrained, and especially CMSSM as 
a solution to the hierarchy problem is strongly disfavored. 

If we motivate new physics models from the reported anomaly of muon g-2, our naive 
expectation is the existence of EW scale new particles because the size of anomaly of muon 
g-2 is as large as one induced by EW gauge bosons in the SM. In this paper, we analyzed two 
examples, where muon has new Yukawa type interactions with new particles. One example 
was a model where right-handed muon has new Yukawa interaction with SU(2) singlet scalar 
and fermion. In this case, in order to explain the anomaly of muon g-2 as well as to satisfy 
the precision EW measurements, we found that the relatively large new Yukawa coupling and 
EW scale new particles are strongly favored. We also noted that the multi-charged particles 
(whose QED charges are large) are favored. 

Another example we analyzed was a model where both right- and left-handed muons have 
new Yukawa interactions with SU(2) doublet scalar and singlet scalar and fermion, so that the 
chirality flip of muon can occur in the internal fermion line of 1-loop muon g-2 diagrams. In 
this case, new contributions to muon g-2 are enhanced, and hence even smaller new Yukawa 
couplings can explain the anomaly of muon g-2. We also showed that the EW observables 
are also well-fitted by the new physics contributions. In this case, the EW scale new particles 
are favored by the data. 

Although these two examples may be only a part of complete models, the muon g-2 as 
well as the EW observables can constrain the QED charge of new particles as well as the scale 
of new particles and the size of new Yukawa couplings. Since the EW scale new particles 
can explain the anomaly of muon g-2 and they can be consistent with the precision EW 
measurements, it is important to study the possible phenomenology at the LHC. 

Since some of new particles have QED charges as well as couplings with Higgs boson, 
the branching ratio of /;.—;■ 77 can be effected. Especially, as we showed, multi-charged 
particles favored by the muon g-2 can significantly contribute to /i — )■ 77, and the effects can 
be sizable. Therefore, Higgs physics will have an impact on models we discussed here. 

We also computed direct production cross sections of new particles motivated by the muon 
g-2. Since models discussed here may not be complete, we could not discuss the complete 
signatures for the models. However, since the number of signal events are not negligible, 
there may be possibility to reveal new physics models for muon g-2 at the LHC. Since we 
do not need the colored particles to explain the anomaly of muon g-2, the search for new 
particles which only have EW interactions will be very important. 

Finally, we comment on the possible completions of models discussed in this paper. The 
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structure of the model where both right- and left-handed muons have new Yukawa type 
interactions is quite similar to that of the minimal supersymmetric standard model. When 
Q2 = —1, the scalar 02 in SU(2) doublet scalar $ and SU(2) singlet scalar cj) correspond 
to left- and right-handed smuons, respectively. The SU(2) singlet fermion x is a Bino, for 
example. A part of a model discussed in Ref. [21], which is motivated by neutrino mass and 
dark matter is also similar to the model we discussed here. Therefore, our approach can 
capture crucial features of new physics models for the anomaly of muon g-2 and hopefully 
create an interesting way to find the physics beyond the standard model. 
Note added 

After we completed this work, we learned that ATLAS and CMS collaboration at LHC 
observed a new particle consistent with the SM Higgs boson at a mass of about 125 GeV |22] . 
Since we have taken the Higgs boson mass to be 125 GeV in our analysis, our conclusion 
does not change. 

Acknowledgments 

This work is supported in part by Grant-in-Aid for Scientific research from the Ministry of 
Education, Science, Sports, and Culture (MEXT), Japan, No. 22540273 (KT). This work is 
also supported in part by Grant-in-Aid for Nagoya University Global COE Program, "Quest 
for Fundamental Principles in the Universe: from Particles to the Solar System and the 
Cosmos" , from the MEXT, Japan. 



26 



A Useful function 

In this paper, Passarino-Veltman functions |9j are defined by 



2,.2e 



Bo{A,B;p) = 167r> 



i(27r)" P - + ie' 



2,.2e 



p^Bi{A,B-p) = 167r> 



i(27r)" - + ie] [(fc + p)^ - m% + ie] 



i(27r)" - + k] [(fc + pf - m| + ze] 

pVi?2i(Ai?;p) + g^'B^^iA^B-p) 



Co(A,5,C;pi,p2) 



2 ,.2e 



IGtt u 

' i(27r)" - + ze][(A; + - m| + ze][(A; +pi +^2)^ - "^c + 
(rfC^ii+P^C^i2)(Ai?,C;pi,p2) 

IQn a / 

y i(27r)" — rn\ + ze] [(A; + pi)^ — m| + ze] [(A; + pi + P2Y — m^c + 

+ P^P^C722 + + KP2 )C^23 + g'^C^,} (A, 5, pi, P2) 



2 ,,2e 



IGvr /i 



i(27r)" [A;2 — + ie] [{k + pi)^ — m\ + ze] [(A: + pi + P2Y — + ze] ' 

(A.3) 



wfiere we use dimensional regularization in 4 — 2e dimensions. 

For convenience, we list explicit forms of some of above functions. 



2x 2 ' 1 - - 



A{m') = m^(- + l-log— ), (A.4) 



B,{A.B;p) = + dxxhg^ ^ \2 ' (A-6) 

B2i{A,B;p) = 77-r-/ rfxx^log^^ ^ 2 ' (A-^) 

3A Jq 11^ 



B22{A,B;p) = ^(ml + ml-^)(^ + l 



1 /"^ 7 r 2/ 2 2 / Nil ml{l — x) + mix — p'^x{l — x) — ie 

- / (ix |m;^(l — x) + m2X — p x(l — x)| log , (A. 8) 



where = 7 ~ 7 + log 47r. 
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